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ABSTRACT 

One of the promising methods to search for life on extra-solar planets (exoplanets) is to 
detect its signature in the chemical disequilibrium of exoplanet atmospheres. Spectra at the 
modest resolutions needed to search for methane, oxygen, carbon dioxide, or water will demand 
large collecting areas and large diameters to capture and isolate the light from planets in the 
habitable zones around the stars. Single telescopes with coronagraphs to isolate the light from 
the planet will have to be 8 m or more in diameter to generate sample sizes with a reasonable 
probability of finding at least one life-bearing planet; interferometers of smaller telescopes can 
overcome some of the limitations but will still need large similarly large collecting areas. Even 
larger telescopes will be needed to detect atmospheric signatures in transiting planets. In all 
cases, the sample sizes increase as the third power of telescope diameter. Direct observation 
using coronagraphs or interferometers are most sensitive to planets around stars with masses 
similar to that of the Sun, whereas transit observations favor low-mass stars near the nuclear 
burning limit. If the technical difficulties of constructing very large space telescopes can be 
overcome, they will be able to observe planets near hundreds to thousands of stars with adequate 
resolution and sensitivity to look for the signatures of life. 
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The detection of more than 200 planets out- 
side the Solar System is a powerful incentive to 
search for extra-terrestrial life. Although extra- 
terrestrial life could take on many guises, econ- 
omy of hypothesis (and a practical approach) im- 
plies that we should first search for signs of life 
similar to those seen on Earth. An advantage of 
using the Earth as a proxy for analysis is that it 
bounds the problem and provides concrete exam- 
ples of signatures subject to passive detection, i.e. 
not requiring signals broadcast by sentient beings. 

In the Earth's present atmosphere, the chem- 



ical components have been altered by life. Ev- 
idence for life on Earth could be detected from 
afar in the spectral signatures of these molecules: 
oxyg en, carbon diox ide, methane, and water va- 
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por (jSagan et al. 19931 ). The rise of photoplank- 
ton and plants on Earth created an atmosphere 
with a large reservoir of oxygen today that requires 
steady production by photosynthesis to maintain 
its present level. If we could study the same 
spectral signatures in the atmospheres of exoplan- 
ets, we could search for signs of life similar to 
some of the earliest and most robust forms on 
our own planet. But the Earth's atmosphere has 
been altered by life several times over the last 
4 Gyr, and there are potentially many signatures 
that co uld indicate the presence of life on exo- 
planets (iTraub fc Jucks 2002t ISeager et al. 20021: 
Kaltenegger et al. 2007[ ). Chemical signatures of 
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life on other planets would revolutionize our think- 
ing about Earth's uniqueness and provide tantaliz- 
ing evidence that we are not alone in the universe. 

Observing exoplanets directly is difficult ow- 
ing to their proximity to the much brighter 
stars that keep them warm. Although tech- 
nically challenging, this problem is well under- 
stood, and there are a variety of strategies that 
can reduce the brightness of the starlight with- 
out diminishing the light from th e planet for 
direct detection (jCuvon et al. 20061: ICash 20061 ) 
or use the star itself as a background source to 
probe the atmosphere when the exoplanet tran- 
sits the face of the star ( Charbonneau et al. 200 j : 
Ehrenreich et al. 20061 ) . The first technique must 
overcome diffraction in the pupil of the telescope, 
a well understood phenomenon, and it is possible 
to characterize the detection problem in general 
terms to understand the kinds of instruments that 
will be needed to study exoplanets and search for 
signs of life. The second technique depends only 
on the photometric accuracy of an observation 
and is easy to calculate for any star. 

Any planet supporting life as on Earth must 
satisfy two broad criteria: (1) it must have surface 
temperatures in the range 273 to 373 K, where wa- 
ter is in the liquid phase, and (2) it must have an 
atmosphere. The first criterion is met if the planet 
is in the habitable zone (HZ) around the star, a 
range of orbital distances where the equilibrium 
temperature for a rotating body is between the 
freezing and boiling points of water. The second 
criterion is met if the planet is rocky and can retain 
an atmosphere; current estimates specify a mass 
between 0.5 to 10 Earth masses. Smaller plan- 
ets will not retain their atmospheres, and larger 
planets accrete gas and become gas giants. These 
criteria are necessary but not sufficient to create 
Earth-like life. Although they are probably far 
too restrictive to encompass all the possibilities 
for other life forms in the universe — or even on 
Earth itself — they are the only ones amenable to 
remote observation with technology that we can 
foresee at present and thus provide a good basis 
for a targeted search. 

There have been many calculations aimed at 
refining our ideas of the habitable zone (e.g. 
Kasting et al. 199^. suitable sa, m ples of stars to 



search dTurnbull fc Tarter 20031 : iTurnbull 200l 
Turnbull et al. 2006h . and the impact of specific 



telescopes on such a search (jAgol 20071 ). and there 
is some disagreement about the likelihood of suc- 
cess depending on the different assumptions used. 
Most authors concentrate on photometric detec- 
tion alone, ignoring the means to search for life 
(e.g. Agol 2007). 

But the utility of photometric searches alone to 
identify exoplanets for subsequent study could be 
obviated by the difficulty of measuring the orbits 
accurately enough to recover them at a later time 
([Brown 20051: iBrown et al. 20071 ). For apertures 
sizes under discussion for the Terrestrial Planet 
Finder (TPF) mission, habitable zones of most 
potential target stars are blocked by the coron- 
agraph. The planets of interest will spend the 
majority of time behind the central obscuration 
of the imaging instrument. To ensure efficient re- 
covery of a newly discovered planet at future ob- 
serving epochs, it will be necessary to estimate 
the orbit to high accuracy from a small number 
of astrometric measurements. This requirement 
implies a lower limit to the aperture size based 
on operational requirements (Brown et al. 2007). 
Achieving adequate astrometric precision may de- 
mand apertures larger than any so far discussed 
for TPF (Brown 2007, personal communication). 
If true, discovery and immediate spectroscopy of 
candidate sources near the stars may be the most 
efficient means of identifying those most interest- 
ing for follow up observations to determine if life is, 
indeed, present, making it essential to understand 
the spectroscopic requirements at the outset. 

The purpose of this article is to derive the main 
scaling parameters for the study of life-bearing 
exoplanets in known samples of stars to under- 
stand the size of the telescopes needed for a ro- 
bust search. We adopt simple but optimistic as- 
sumptions to bound the problem and find the min- 
imum size for survey telescopes. Using only the 
lowest order approximations and assuming "best 
case" observing conditions allows robust conclu- 
sions about the scale of facilities needed to tackle 
the search for life. 

The main premise is that direct photometric 
detection of exoplanets in a band where the exo- 
planet atmosphere is free of chemical signatures 
cannot be the endpoint of any mission to search 
for life-bearing planets; spectra of the atmospheres 
will be the major advance of observing the planets 
directly. Moreover, the rapid increase in the num- 
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ber of exoplanets discovered to date suggests that 
finding terrestrial planets will be easiest with in- 
direct methods, such as observing radial velocity, 
photometric or astrometric variations in the host 
stars, and the real thrust of direct observations 
will be to search for signs of life. 

2. Detecting Earth 

The Earth's flux density viewed from a distance 
is reasonably well constrained by its size and the 
requirement that its equilibrium temperature is 
determined by external illumination at the tem- 
perature of the Sun. The same statements ap- 
ply to exoplanets in the habitable zones around 
distant stars. An important result is that the 
flux densities of exoplanets in the habitable zones 
across a broad spectrum are roughly independent 
of the stellar luminosity and temperature. 

If the exoplanet is a spherical body with radius, 
Tp, and temperature, Tp, and it is viewed at a 
distance, D ^ Vp, it's thermal emission is at most 
that of a blackbody, B: 



(1) 



The total emission is ^irrpCrTp. 

The reflected light is in principle a rather com- 
plicated function of the planet's orbital phase 
and the det ailed scattering pro perties of its atmo- 
sphere (e.g. Seager et al. 200ci[) . For the purposes 
of this paper, it will be sufficient to assume the 
planet scatters like a Lambert sphere viewed at 
quadrature, the viewing angle a = 90°, with a 
fraction 1 — a{v) of the light absorbed in the pro- 
cess. If it is illuminated by a star of radius, r*, 
with an effective temperature, T*, at orbital dis- 
tance, from the star, the observed flux density 
of the reflected light is: 

FA^) - --^<^(a)aH-^^B(^,T,), (2) 

= la{v)^^7:B{v,T^), (3) 

where the Lambert sphere has a geometrical cross 



section of | 



± sm el- 



and a phase function q!)( a) 

(TT — Q) cosa) = i for a = ^ ( Russell 1916t 
Seager et al. 200Clf) . The unabsorbed fraction, 
a{v), will in general depend on frequency. For 



the Earth, it varies between 0.7 near 0.3 ^m 
to less than 0.3 longward of about 0.6 ^m and 
depends on variable factors such as cloud cover 
and exposed land i nass ([Kasting et al. 19931: 
Turnbull et al. 20061: iMcCullough 2006h . 
adopt a{v) = a = 0.5 for the rest of this paper 
to produce a bright planet that should be most 
easily observed. As shown below, uncertainties in 
a{i') will not affect the estimates in this paper sig- 
nificantly. The actual spectra would also include 
radiative transfer effects in the atmosphere and 
show absorption lines from molecules (see Fig. [1]). 
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Fig. 1. — The lower line is a synthetic spectrum of 
Earth viewed from 10 pc to illustrate the spectral 
features in the Earth's present-day atmosphere. 
The upper line is a combination of the solar ir- 
radiance and a Planck function at 286 K. 

Combining ([1]) and ([3]) with the assumptions 
above, the flux density of the planet is: 



Fp{i^) 



p 

i?2 



1 K 
3Rp^ 



B{v,T,) + B{v,Tp)\. (4) 



This relation neglects some interesting details, 
such as the variable atmospheric transmission and 
radiative transfer effects as well as adopting as- 
sumptions about the albedo and phase function 
that are at best educated guesses. It turns out 
that these simplifications do not affect the overall 
energy balance enough to change the conclusions 
of this article, but they are, nevertheless, crucial 
for understanding the chemical composition of the 
planet's atmosphere. 

At 10 pc, the apparent magnitude of the Earth 
at visual and near infrared wavelengths from ^ 
is about 29(AB), ~ lOnJy, too faint for spectra 
with any existing telescopes. 



3 



The principal sources of noise in detecting exo- 
planets will be local zodiacal light, any equivalent 
exo-zodiacal light in the other planetary system, 
any residual light from the star not cancelled by 
the instrument, and the planet itself. Since we 
are interested in the limits of nature, we assume 
that starlight is completely eliminated by the tele- 
scope and instrument and consider only the other 
contributions. 

The zodiacal light has the same spectrum as 
the Sun at wavelengths shorter than about 3 /im, 
but reduced by the optical depth of the zodia- 
cal cloud in the direction of observation, and it 
is approximately a blackbody function at thermal 
infrared wavelengths with a temperature that de- 
pends somewhat on the viewing angle with respect 
to the Sun. The important point is that both the 
zodiacal light and exo-zodiacal light light will be 
uniform surface brightnesses to a good approxi- 
mation whose observable fluxes depend only on 
the area-solid angle product (etendue) of the tele- 
scope and not the distance to the source. The 
spectrum of the zodiacal light combines a visual 
optical depth, ^ 10^^, and thermal emissivity, 
^z{v) ^ 0.5, both strong functions of ecliptic line 
of sight, to obtain the brightness of the zodiacal 
light as: 



t,[B{v,Tq) + e,B{i^,T,ody)]. (5) 



It is straightforward to calculate the exact zodiacal 
emission in any direction from models in the liter- 
ature (e.g. Leinert et al. 1998), and as a practical 
matter we will adopt a viewing angle at solar elon- 
gation 90° and ecliptic latitude 45° as "typical" of 
a large survey seeking to minimize zodiacal back- 
ground: /zody('^) is about 0.3MJysr~^ at 0.5 /xm 
and lOMJysr"^ at lOyitm. Figure [2] compares the 
photon rates seen by 4 and 8 m telescopes from 
an Earth 10 pc away with the zodiacal light for an 
assumed resolving power of 100. 

It is impossible to know how strong the contri- 
bution from the exo-zodiacal light is. It is easy 
to find stars wi th very large quantiti es of inter- 
planetary dust ( Beichman et al. 2006f l and stars 
without detectable dust emission, so this term is 
one of the larger uncertainties in our estimates 
as stressed by Agol (2007). The best case is if 
there is no exo-zodiacal light at all. Alternatively, 
if the other planetary system is identical to the 
Solar System, the amount of exo-zodiacal light is 



just twice the local value. The easy way to un- 
derstand this factor of 2 is to realize that we view 
the exo-zodiacal cloud through its entire thickness, 
whereas we reside at the mid-plane of the local zo- 
diacal light (within 3°) and thus look out through 
only half the optical depth in any direction. The 
exo-zodiacal cloud will also be well resolved by any 
telescope that can block the starlight and still ob- 
serve the planet. For the moment, we assume the 
contribution of the exo-zodiacal light is negligible. 
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Fig. 2. — The photon rates of an exo-Earth viewed 
from 10 pc distance with 4 m and 8 m telescopes, 
together with the photon rates for the local zodi- 
acal light for a diffraction-limited telescope. 

The noise of an observation will be dominated 
by the (Poisson) statistics of the incoherent pho- 
ton stream from the planet and the zodiacal back- 
ground. Denoting photon rates as Nj{p) and 
N^{z), respectively, the signal to noise ratio, SN, 
after observing time, t, is: 



SN = 



1/2 ' 



1 Aiy At^iFpjiy) 



(6) 



h V ^^tel (^;,(i^) -f f^tcl4(!^))' 

(7) 



with the area and solid angle of the telescope 
denoted and iltch a-nd the overall detec- 
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tion efRciency rj. For a circular telescope that is 
difFraction-limited, the etendue, Atciritoi, is just 
equal to the square of the observing wavelength, 
A^. In the limit where the planet is far away from 
the observer (us), the photon noise is dominated 
by the zodiacal light regardless of whether it is 
local or around the exo-planet, since its contribu- 
tion is independent of both source distance and 
telescope diameter, but the photon flux from the 
exo-planet is proportional to D^^. Equation ([7]) 
may then be written as: 



SN 



dtci 

IT 



(8) 
(9) 



The expression in ^ shows the explicit depen- 
dence of SN and observation time on telescope 
diameter and source distance in the faint limit. If 
exo-zodiacal light is comparable in brightness to 
the local value, SN diminished by a factor of or- 
der VS- It is useful to note that even though the 
amount of exo-zodiacal light is uncertain, it intro- 
duces only a modest uncertainty into the signal-to- 
noise calculation considering many of the other un- 
certainties inherent in the planet detection prob- 
lem. However, we will assume it is zero in sub- 
sequent calculations to consider the best case sce- 
narios for taking spectra of planets. 

Recasting (jS]) in terms of observation time for 
Earth at a distance is especially useful for the 
planet detection problem. Using ^ and ([5]), Fp 
and Iz can be written in terms of assumptions 
about the star, planet, and observing system pa- 
rameters. To simplify later calculations, we define 
SNw = SN/10, Two = v/Aiy/100, dg = dtci/8m, 
Dio = D/lOpc, and ^24 = t/24hr: 



SN^oV'^Two 



assuming an observation wavelength A = 1 /im, 
justified below. Equation pH)) shows that the time 
for an observation increases as the fourth power 
of distance over telescope diameter when observ- 
ing a point source against a uniform brightness 
background in the photon limit. Figure [3] plots 
the times need to take low signal-to-noise spectra 
of an Earth at 10 pc with 4 and 8ni telescopes. 



This strong dependence means that observations 
of Earth-like exoplanets with a telescope of a given 
size will be effectively limited to some maximum 
distance: 



= 40 pc ds 
« 30 pc c?8 



77*24 



T^woSN?Q 



1/4 



(11) 
(12) 



where we adopted r/ — 0.25 and set the other pa- 
rameters to unity for the purposes of later sim- 
plification. This equation applies to observations 
of refiected light at l^m from the planet and is 
strictly valid only for distances where the planet is 
fainter than the zodiacal light. For bright objects, 
the source flux itself also adds to the noise further 
increasing the required integration time. The lim- 
iting distance for longer wavelengths using ^ will 
have the same functional dependence but a smaller 
distance normalization given by the relative in- 
crease in the integration time at the longer wave- 
length (Fig. [3]) to the negative one-fourth power. 

One of the interesting results of this calcula- 
tion is that the observation time near the peak 
of the thermal radiation, 10 /im, is about the 
same as at shorter wavelengths despite the enor- 
mously larger background radiation. The reason 
is because the background noise increases only as 
the square root of the photon rate, and very high 
photon rates reduce the relative noise from the 
background very rapidly, at which point the high 
photon rate from the exoplanet is easy to detect. 
The only real disadvantage of thermal infrared ob- 
servations is the poor diffraction-limited resolution 
of a single filled aperture, making it difficult to 
detect planets at small angular distances from the 
stars. 

The choice of 1 fim for pop is appropriate for 
any observations with a single-telescope, where the 
ability to eliminate starlight is tied to the tele- 
scope diameter as discussed below. Observations 
in the thermal infrared at wavelengths longer than 
a few microns encounter enormous challenges of 
starlight suppression at small angular distances 
from the star, making single telescope observa- 
tions unworkable. However, multiple telescopes 
used as an interferometer or the combination of a 
single telescope with a distant occulter would sep- 
arate the observational sensitivity considered here 
from the suppression of starlight. We will con- 
sider both cases in the following sections to show 
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the independent effects of sensitivity and starlight 
suppression. 
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Fig. 3. — Integration time to achieve S/N=10 
for each spectral element for a practical 4 m and 
8 m telescope at a resolving power of 100 as- 
suming only background from the local zodiacal 
light; integration times would increase for any exo- 
zodiacal emission. The overall detection efficiency, 
77, is 0.25 for this figure. 

The assumptions used thus far assume a per- 
fectly diffraction-limited telescope where the fore- 
ground radiation is dominated by zodiacal light. 
These conditions are easy to meet for a space 
telescope; they are extremely difficult to meet for 
any telescope on the Earth's surface. Even if the 
distortion of the wavefront from the Earth's at- 
mosphere could be entirely compensated by an 
adaptive optics system, the atmosphere radiates 
strongly and absorbs light from the exoplanet in 
all the lines of interest in Fig.[l]making the charac- 
terization problem almost impossible for a ground- 
based telescope. 

The practical effect of pUj) is to limit the use- 
ful distance for observations of an exo-Earth with 
a single telescope. It takes about 24 hr to take a 
spectrum of Earth at a distance of 40 pc with an 
8 m telescope, a distance which shrinks to 27 pc, 
if the detection efficiency is 25%, typical for a 
good spectrograph. For a 4m telescope, the lim- 



iting distances for a 24 hr observation time are 20 
and 13.5 pc, respectively, since the distance scales 
linearly with telescope diameter for this problem. 
These sizes are larger than the current generation 
of space telescopes. It is also useful to note the 
size of telescope below which the zodiacal light 
begins to dominate the light from the planet as 
a function of the distance to the planetary sys- 
tem, Fpiy) < ^teilzodyii^) in the denominator of 
d?]). Using dH) and ([5]) at a wavelength of 1/im, 
c^tei < 3.25m (Z?/10pc). At thermal infrared 
wavelengths, zodiacal light dominates that light 
from any interesting exoplanet for all telescopes 
less than a few tens of meters in diameter. 

An interferometer can be treated as a telescope 
with a sparsely-filled aperture. The total collect- 
ing area is just the sum of the individual collecting 
areas. Each telescope collects background radia- 
tion over a solid angle set by the diffraction-limit of 
each individual aperture, the "main beam," modi- 
fied by the interference pattern created by the spa- 
tial separation of the telescopes. 

One way to think about this problem is to note 
that the solid angle is proportional to the Fourier 
transform of the pupil-plane (u-v) coverage, which 
must be equivalent to the pupil-plane coverage of 
a single telescope with an area equal to the sum of 
the areas of the elements. Therefore, if the sum of 
the collecting areas of the interferometer elements 
is Atot, the sensitivity of a single observation can 
be calculated using ([5]) with dtoi = v^A^ot/TT- If 
a complete aperture synthesis is required to take 
spectra with the equivalent resolution of a single 
telescope with diameter, c?i, the total time will 
be ndl / (AAtot) times the single observation time. 
Full synthesis may not be needed if the position of 
the exoplanet is known from some other observa- 
tion. 

Sample volumes of the local neighborhood in- 
crease as the cube of the distance. Because there 
is a premium on maximizing the number of survey 
stars that can host potentially habitable planets, 
these results show a strong preference for large 
telescopes or interferometers with large collecting 
areas. The preferred size will depend on the den- 
sity of suitable stars in the Solar neighborhood 
and their luminosity distribution to set the scale 
for the size of the habitable zones, the subject of 
the next section. 
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3. Sample Sizes for Direct Observations 

The interesting candidate stars will be those 
whose habitable zones can be studied after the 
light from the star is suppressed to a level be- 
low the background. Because the starlight must 
be essentially eliminated — the Sun is ~ 10^" times 
brighter than its reflected light from the Earth — 
starlight cancellation is incredibly demanding. It 
is generally a far higher technical hurdle than de- 
tection for exo-planet observations, and it is the 
reason that articles on this subject often underes- 
timate the importance of background-limited sen- 
sitivity when considering exo-planet studies. 

For a single telescope, starlight suppression re- 
quires a coronagraph that rejects light over a small 
area in the focal plane, while allowing observation 
of the light from outside this area. The smallest 
area that can be darkened depends on the lim- 
iting resolution of the telescope which can never 
be smaller than the diffraction limit of the pupil. 
An external occulter suppresses starlight by using 
a distant screen decoupled from the telescope to 
create a shadow of the starlight outside of which 
exoplanets may be studied. An interferometer is 
essentially a sparse aperture telescope that can be 
phased to null starlight like a coronagraph over an 
angular area the depends (inversely) on the sepa- 
ration of the elements. 

The inner working angle, 9iwA, is defined 
as the smallest angular separation between a 
bright star and a much fainter planet at which 
the planet can be studied with good cancel- 
lation of the starlight. For a single telescope 
with a circular pupil, the diffraction p attern 
is an Airy function ( Born &: Wolf 19991 ) . and 
OiWA cannot be less than the distance to the 
first null at 1.2A/(itoi- Most realistic corona- 
graphs achieve inner working angles more than 
twice this angle, even w ith nearly optimal de- 
signs ( Guvon et al. 20061 ). The angle of the first 
null is typically many tens of milli-seconds of arc: 
1.2A/(itci — 0.21 arcsec (A/1 ^va) (dtci/l m)~^. 

For a telescope-occultcr combination or an in- 
terferometer, the inner working angle is not con- 
strained by the diameter of the light collecting 
telescope (s). An occulter produces a fixed inner 
working angle determined by its size and distance 
from the telescope. The inner working angle of an 
interferometer depends on the maximum distance 



between individual elements (telescopes) and will 
vary with position angle on the sky, depending its 
configuration. The resolution of an interferome- 
ter varies with wavelength in the same manner as 
a coronagraph, allowing us to characterize the in- 
ner working angle as a fraction of the diffraction- 
limit of the effective diameter of the telescope, 
OiwA = /el-2A/dtci for any system, with /g > 1 
for a coronagraph, /e <C 1 for an interferometer, 
and fg (A) a function of wavelength for an occulter 
but generally less than one for the wavelengths of 
observation. 

In all cases, the size of the star sample available 
for study will increase as the inner working an- 
gle decreases, making the smaller habitable zones 
around numerous low-luminosity stars observable 
at any distance. Put differently, as 9iwA shrinks, 
the minimum mass of main sequence stars avail- 
able for study decreases, and the sample size grows 
from the inclusion of low mass stars. 

We can simplify this problem by deriving scal- 
ing relationships for the size of the habitable zones 
around stars as a function of stellar parameters 
and characterizing the samples in terms of the 
present day mass function of stars that gives their 
number density as a function of stellar mass. The 
relative advantages of the different star suppres- 
sion techniques will then emerge as a comparison 
of sample sizes for specific observing techniques. 

Almost all the interesting spectral features in 
the Earth's atmosphere are longward of 0.7 /zm 
(Fig. [TJ) A spectrum should extend to at least 
1 /im to be analyzed for evidence of disequilib- 
rium chemistry indicative of life, and longer wave- 
lengths bring in even more interesting features. 
We adopt A = 1 /^m as the minimum working 
wavelength to look for signs of life in spectra and 
note the desirability of observations to look for sig- 
natures of ozone, methane, and water that have 
strong bands throughout the thermal infrared be- 
yond 10 jjwi, where integration times are near a 
second minimum as seen in Fig. [S] 

The size of the habitable zone depends on the 
stellar luminosity. The total thermal emission, 
47rr^crTp , must equal the total amount of absorbed 
starlight, (1 — a)i*/(47ri?p). The planet's temper- 
ature, Tp, is bounded between the freezing and 
boiling points of water, 273 K < Tp < 373 K, a 
necessary but not sufficient condition for life as it 
is on Earth. We assume for this calculation that all 
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the starlight is absorbed (a = 0) to maximize the 
habitable zone radii. This condition is sufhcient to 
compute the inner and outer habitable zone radii 
for a star of luminosity, using general formulae 
for the orbital radius, Rp, of a rotating planet at 
equilibrium temperature, Tp-. 



167rcrr4' 



(13) 
(14) 



The inner and outer radii of the habitable zone 
are obtained from (fH|) with Tp ~ 373 and 
273 K, respectively. Effects such as greenhouse 
warming or low absorption of starlight via high 
albedo will change t hese radii by modest factors 
(|Kasting et al. 1993f) m opposite directions; for 
example, th e effective t emperature of the Earth 
is - 255 K (jOrton 2000l ) and the albedo is - 0.5 
but the greenhouse effect keeps the surface warm. 
These corrections are not large enough to change 
the overall conclusions of this paper, however, and 
we will subsequently ignore them in the interest 
of simplicity. 

Using to compute the radii bounding 

the habitable zone, OjwA determines the maxi- 
mum angular area that can be blocked to achieve 
starlight suppression for a star at a known dis- 
tance and luminosity. By examining large sam- 
ples of stars, it is possible to determine how many 
have habitable zones that can be observed with 
a telescope or telescope system (interferometer or 
telescope and occulter). 

To the accuracy required for the treatment 
in this paper, we express the stellar parameters 
in terms of the normalized stellar mass, m = 
m*/M©, as follows: 



Lq m 
Tq m 
To 



(15) 
(16) 
(17) 



The approximations are good to about 25% for the 
local sample; for the Zero Age Main Sequence, the 
coefficients in front of the three equations would 
be 1.2^0, O.96r0, and I.ITq, respectively, and the 
exponents would be 3.8, 0.92, and 0.49, respec- 
tively. 



These relationships allow us to recast the equa- 
tions for observing the atmospheres of exoplanets 
in terms of stellar mass. In particular, we can sep- 
arate the impact on sample size of the collecting 
area from inner working angle, allowing any ob- 
serving system to be optimized when trading off 
these two parameters. 

The star is relevant only to the size of the hab- 
itable zone through , the planet always having 
approximately the same spectrum (we will adopt 
T(j, — 285 K to correspond to earlier sections): 

2 

I I I \ 

R„ 



1.0 AU 



(18) 



(19) 



The size of the habitable zone is proportional to 
■\/L^ oc from (jl3p and p5p . For any inner 
working angle, there is a minimum mass star for 
which the habitable zone might be observed at any 



distance, Rp{m-a 



OjwaD. We will specify the 



inner working angle as a multiple of the diffraction 
limit of a circular pupil, 9iwa — /9(l-2Ap/(itci) 
and for generality define the equivalent diameter 
of the telescope system in terms of the collecting 
area, dtoi = A/4Atci/7r. Then, the maximum dis- 
tance to a star of normalized mass, m, where the 
habitable zone can still be resolved is: 



DiwA{m) 



Rp/9 



■p/OlWA 



SOpcdsf, 



-1 I 



(20) 
(21) 



For any star, the maximum distance for spectro- 
scopic observations will be the lesser of DjwA or 
Ds, the limit from the integration time in (|12[) . 
The minimum stellar mass whose habitable zone 
is observable at distance, D, is: 

1/2 



fe D 
ds 30 pc J T( 



(22) 



The present day mass function, ^(m), of local 
stars above the nuclear burning limit is well ap- 
proximated by tw o power law^j in stellar mass 
(|Reid et al. 2002f) : 



0.02 pc" 



(23) 



^ Other choices are possible l|KrouDa 200lh . but this choice 
fits both the 8 pc sample and the more extended sample of 
solar mass stars to 25 pc. 
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with q;(0.08 < m < 1.0) = 1.1 w 1 and a{m > 1.0) 
= 5.2 « 5. The PDMF is defined such that the 
number density of stars in the mass range m to 
m + dm is ^{m)dm. The large exponent for the 
upper end of the mass range describes the actual 
number density of stars above a few solar masses 
as opposed to the initial mass function which is 
derived by correcting for stellar evolution. 

The sample size distribution for direct observa- 
tions follows from (j23p using the maximum volume 
for each stellar mass: 

dNd/dm = (24) 

where D is either Djwa or Dg- The stellar mass 
for which Dg becomes the limiting distance de- 
pends on the rejection factor, /g, and the place 
in the habitable zone where the planet resides, 
governed by Tp. If we assume a perfect corona- 
graph observing Earth, then fg{Tp/T,^) « 1, and 
Ds < Di\YA for m > 1, exactly where the expo- 
nent in ^(to) changes value. The total sample for 
all masses above 0.1 M© is: 

Nd = ^{ / D^jwji{m)^{m)dm 

+Dl ^{m)dm^ (25) 

«^ 1900 dl (26) 

The sample size depends strongly on the effec- 
tive diameter of the telescope, inner working angle, 
and working wavelength in the low-mass limit. It 
depends weakly on the overall detection efficiency, 
integration time, signal-to-noise ratio, and spec- 
tral resolving power. 

3.1. Samples for Coronagraphs 

Direct observations of planets with corona- 
graphs tend to be limited by the inner working 
angle and benefit from the largest possible hab- 
itable zone. Evaluating l|26p with fg = 2, the 
sample size as a function of telescope diameter is: 

Nd « 280 (^^) dl (27) 

An 8 m telescope with an excellent coronagraph 
should have almost than 300 stars around which 



direct characterizations of Earth-like planets could 
be made spectroscopically, assuming every star 
where suitable for observation (e.g. no binaries). 
But the dependence on the third power of tele- 
scope diameter shows that the best one could do 
with a 4m telescope is about 35 stars, and this 
is the total sample without correction for multiple 
star systems and stars that might not be suitable 
for life. On the other hand, a very large telescope, 
16 m, say, could study thousands of stars to char- 
acterize terrestrial planets. 

3.2. Samples for Compound Telescopes 

The samples can grow rapidly when fg is less 
than 1 . Use of an external occulter should produce 
fe ^ 1, and an interferometer could give /e <C 1, 
although the instantaneous coverage of the {u, v) 
plane is sparse and would be of limited value unless 
the precise location of the planet were well known 
at the time of the observation. 

Evaluating ([25]) for the case fg = 0.1 and Tp = 
T® yields: 

Nd ~ 3300 d|, for /e= 0.1 (28) 

This is the sample a good interferometer might 
access. An occulter can reduce fg to less than 1 
(e.g. Cash 2006), but the sample size grows rather 
slowly with decreasing fg owing to the sensitivity 
limit. The typical size of an occulter to produce 
fg<l coupled with its distance from the telescope 
to produce a small inner working angle will make 
a practical occulter difficult to use. 

These examples show that the gain from us- 
ing a larger telescope comes about mainly from 
the rapid increase in sensitivity with aperture size, 
d^j, and the rapid increase in sample size owing 
to the larger volume probed: TV ~ ~ d^^j. 
The use of an external occulter helps mitigate 
the technical demands of creating a nearly per- 
fect coronagraph, although it will drive some very 
strong technical demands of its own — alignment 
of a ~ 60 m structure at 400,000 km from a space 
telescope and moving it around to sample different 
stars will be no mean feat. Large telescopes will 
still have an enormous advantage for searching for 
signs of life on exoplanets. 
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4. Spectra of Transiting Exoplanets 

It is possible to distinguish very slight differ- 
ences in the light curve at different wavelengths 
as an exoplanet transits the face of its central 
star. These differences result from transmission 
spectrum of the star through the planet's atmo- 
sphere making the apparent size of the planet 
larger at wavelengths where the atmosphere is 
opaque. Charbonneau et al. (2002) and Vidal- 
Madjur et al. (2003, 2004) demonstrated the util- 
ity of this method by detecting sodium, hydrogen, 
oxygen, and carbon in the extended atmosphere of 
the transiting exoplanet, HD 209458b, using the 
STIS spectrograph on 2.4 m Hubble Space Tele- 
scope. Tinetti et al. (2007) used the Spitzer 
Space Telescope in a similar way to detect H2O 
in HD 189733b. Transits provide an alternative 
way to look for the signa tures of life in exoplanets 
(|Seager fc Sasselov 2000l ). 

As a planet passes across the face of a star, 
the starlight dims by approximately the ratio of 
the area of the planet to that of the star, Ap « 
7rrp/7rrJ, where we ignore the modest effect of limb 
darkening. If the scale height of the atmosphere is 
hatm and an atmospheric absorption feature has 
the effect of blocking all light to x scale heights 
above the surface, there will be a small difference 
in the depth of the light curve viewed at the wave- 
length of the feature: 

A„,^-A, = + (29) 



2xh. 



atm'^p 



(30) 



where the approximation assumes xhatm ^p. 
For the Earth crossing the face of the Sun 
seen from afar, this an excellent approximation: 
hatm ~ 9 km, Tp w 6400 km, and Aatm - Ap = 
2 X 10^^ for X = 1. The scale height of the atmo- 
sphere depends on the density and temperature 
of the planet through the usual thermodynamic 
relation: 



(31) 
(32) 



where g is the gravitational acceleration at the 
planet's surface, p. is the mean molecular mass. 



and p is the mean density of the plane10. Recast- 
ing ([50]) in terms of the density: 



27rG'^prJ 



(33) 



This result shows somewhat non-intuitively that 
the atmospheric signature for a planet in the hab- 
itable zone is independent of the planet's size; to 
first order, it depends only on the density and tem- 
perature, the latter of which lies between the freez- 
ing and boiling points of water, thus bounding the 
range of possibilities. The stronger gravitational 
field of larger planets decreases the scale height to 
compensate for the increase in the planet's radius. 

The number of scale heights for the effective 
cross section, x, depends on the strength of the 
absorption lines. In the Earth's atmosphere, the 
strength of the absorption lines is normally ex- 
pressed in terms of the vertical optical depth from 
the ground to space. For a molecule whose number 
density, n^, as a function of height, z, is exponen- 
tial with scale height, hatm, and absorption cross 
section per molecule, ct^, the vertical optical depth 
is '^0 = niaiexp{-z/hatm) dz = riiaihatm- 
Then the optical depth through the atmosphere 
along a line with an impact parameter x scale 
heights above the surface is given by: 



t{x) 



To 



2roe< {r'p + x)Ki{r'p + x), 
ToV2^(r'+xV 



dy (34) 

(35) 
(36) 



where r'^ — rp/hatm, Ki is a modified Bessel func- 
tion, and the approximation assumes -I- a; ^ 1. 
The value of x that gives the effective atmospheric 
cross section for any absorption line seen in tran- 
sit is the solution to (|36p that gives t{x) = 1. For 
r'p ^ X, the solution is: 



0.92 
4.20 



1 



■In 



2 \ hatm 

In(ro) for r'p 



In(ro) (37) 
6400 



(38) 



^The density of a planet increases slightly with mass a s the 
interior equation of state changes l lSeager et al. 2007t) . but 
the effect is small enough to be ignored at this level of 
approximation. 
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The lines shortward of 1 fim in Fig. [T] have 
vertical opt ical depths, tq, between 0.1 and 1 
meaning x is between 2 and 4. There 
are very strong H2O lines at wavelengths of a few 
microns, in principle giving x ~ 10, although for 
the Earth itself, water vapor is confined to the 
tropopause within a single scale height, so the lack 
of vertical mixing means that x(H2 0) ~ 1 despite 
the large tq. The important result of (|37|) is that 
the effective cross section within an absorption line 
depends on the logarithm of the line optical depth; 
even order of magnitude uncertainties in the total 
atmospheric density affect the transit depths mod- 
estly. 

Detecting the absorption lines in a transit ob- 
servation means achieving very high photometric 
accuracy in the difference spectrum in and out of 
the line. The best photometric accuracy that can 
be achieved during a transit is if all the starlight is 
detected with perfect efhciency and the only noise 
source is the statistical fluctuation of the photon 
rate as assumed in ([7]). Using the formalism above 
for measurement time, tobs, and assuming the dif- 
ference of two measurements: 

/ o \ 
cn^ I 2 \ 



^Nytobs 

8h~ 



-1/2 



(39) 



(40) 



The total observation time will be the time of a 
single transit times the number of transits. Viewed 
from a great distance, the maximum transit time 
is just ttrans = ^.T^jvp = 2r^^ R^jGM^, and the 
time for an average transit (i.e. not across the 
equator) is just one half this value. The number 
of transits observed during the lifetime of a space 
mission, t^, is just np = tm/P, with np truncated 
to an integer . The total observation time is: 



tob 



nptf 



_ 1 

— t-m ^ 
TT Hn 



Rr, 



tn 



2 fT, 



(41) 
(42) 

(43) 

(44) 



We used to express Rp in terms of the temper- 
atures to yield another non-intuitive result: the to- 



tal observation time decreases as the stellar effec- 
tive temperature increases. This effect combined 
with the lower contrast for a transit is more than 
enough to offset the advantage of having a brighter 
star with a higher photon rate: 



A, 



3kT^dtei 
47rG/Iprjr* 

ikT^dtei 



(45) 



1 Aiy 

h V 

1 



1/2 



1/2 



(46) 



-tr. 



nr; 



2/W/3 



1/2 



c2 exp(^)-l 



(47) 



where we approximated the flux density from the 
star as a blackbody in (jiTj) . 

Using the same approach as in ^j3] to express 
the stellar quantities, and r^, in terms of the 
stellar mass, setting x = 4 from (|38p . adopting 
jl = 30amu and p — 5gcm^'^, the signal-to-noise 
ratio for a transit observation is: 



SNfr = 240 5^ 



-1.5 



X (exp(2.5m-"-5) - 1) (48) 

The signal-to-noise ratio for transit observa- 
tions of terrestrial planets in the habitable zones 
around stars increases rapidly with decreasing 
stellar mass: small stars are strongly favored over 
massive stars at all distances, both because the 
transit depths are larger and the observation times 
are longer, more than offsetting the increased 
photon rate with increasing stellar temperature. 
This is a general result that also applies to the 
photometric detection planetary transits, not just 
the highly challenging observation of atmospheric 
lines. Notice also that the signal-to-noise ratio de- 
pends linearly on ds/D, a slower rate than for 
background limited photometric observations of 

The challenge of using transits to detect atmo- 
spheric features is illustrated by evaluating (|40|) 
for the Earth/Sun system at a distance of 10 pc. 
A single transit of the Earth across the diame- 
ter of the Sun takes 16 hours during which an 
8 m telescope would capture 3 x 10^^ photons at 
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a spectral resolution of l/F — Aiy/iy — 0.01 at 
1 fim. The fractional uncertainty in the differ- 
ence in light curves at two frequencies would be 
o'N-j/N^ — 2.6 X 10^^, yielding a signal-to-noise 
ratio of 3 on for a feature with T{hatm) ~ 1 
and X — such as the oxygen lines. Since this 
signal-to-noise ratio is too small for most features 
of interest and the required telescope is large, it 
will probably be impossible to detect signatures of 
more distant Earths using transits at these wave- 
lengths. However, large telescopes might be able 
to detect deep water features with T{hatm) ~ 10"* 
at wavelengths of a few microns with sufficient col- 
lecting area, if the water vapor extends throughout 
the atmosphere, unlike the situation on Earth. 

The greater complication with this technique, 
however, is that the orbital plane of the planet 
must be aligned with the line of sight for a transit 
to be observed. The a priori probability of align- 
ment, assuming orientations uniformly distributed 
within the celestial sphere and requiring only that 
the planet graze some part of the stellar surface, 
is: 



maximum near m — 0.12 and falls off very steeply 
above m=l. 

One other quantity of interest is the expected 
number of stars in the sample in which the orbital 
planes will align so that the planets will transit 
the face of the star. Using (|50p with the formalism 
above, the expected number is: 



dNet I dm 



f{m)dNt/dm 

2di (^'' 



(exp(2.5m-o-5) - 1 



,3/2 



(52) 



(53) 



with as — 6.5 for m < 1, or as = 10.5 for m > 
1. The number of observed transits is so heavily 
weighted toward low mass stars that the technique 
favors those near the nuclear burning limit. 

Integrating ([5T|) and ((53)) from m = 0.1 to oo 
yields the total sample size and expected number 
of in the sample in which the orbital plane of the 
planet is aligned for transits: 



/ 



Rp ' 



(49) 
(50) 



where we have used p4)) to derive ([50|) . The frac- 
tion, /, is much less than unity even for cool stars, 
meaning we would need many candidate stars to 
ensure that even a few transiting planets would be 
observed. 

4.1. Samples for transit observations 

The maximum distance for transit observations 
of a star is derived from (|48| by solving for Df. 
Then, using (p3)) for the PDMF, the sample size 
distribution, dNt/dm, is just (47r/3)£'f (m)^(r7i): 



dNt/dm 



1000 ( ^ 

m 



(exp(2.5m 



-o.5^ 



1 



,3/2 



(51) 



with a2 = 5.5 for m < 1, or a2 — 9.5 for m > 
1. The transit technique strongly favors the lower 
mass stars and hot planets near the inner edge of 
the habitable zone. Equation ([5T|l has a strong 



Nt 

Net 



500 d^f^ 



(54) 
(55) 



The final two columns of Table [3] contains the to- 
tal number of candidates and expected number of 
observed transits from (|54p . and ([55]) . 

Notice that the assumptions used to arrive at 
()55p are for the best case of maximum transit time 
and perfect detection efficiency, and there is even 
a slight increase in the observing time if the num- 
ber of transits is not truncated to an integer; this 
result is the best we can do. 

Figure [4] shows how the samples for direct and 
transit observations are distributed across spectral 
type. These are complementary techniques, with 
direct observations favoring stars of about 1 M© , 
and transit observations favoring low mass stars 
near the nuclear burning limit. It will be useful 
to use both techniques to characterize the atmo- 
spheres of planets around stars of all masses. 

These results again point to the great advan- 
tage that large telescopes bring to the study of 
transiting exoplanets. An 8 m telescope should 
see of order 5 transiting planets in the terrestrial 
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zones around low-mass stars in the Solar neigh- 
borhood, if planets are common, and a 16 m would 
yield about 40. The actual number of candidates 
needed to ensure the detection of a life-bearing 
planet depends on how commonly they occur, a 
topic treated in the final section. 

Sample Distributions 



M6 M5 M4 Transits 
M7 



8m telescope 



M4 Direct 




FO 



0.3 0.5 

m,/Msun 



Fig. 4. — The distribution of sample stars vs. 
spectral type that could be observed with an 8 m 
telescope direct detection with coronagraphs from 
([M)) with fe — "2 (blue) and through transit obser- 
vations from ([nj) (red). 



5. Samples from Star Catalogs 

There have been several approaches to deriving 
samples of stars suitable for direct detection of ter- 
restrial planets. One approach is to take all stars 
within a complete nearby sam ple, such as those 
within 8 pc of the Solar System (jReid et al. 20041 ) 
and exclude only the multiple star systems and 
white dwarfs that would be difficult to observe. 
Another approach is to apply criteria that maxi- 
mize the likelihood that some evolved life forms 
might exist by excluding in addition the very 
young stars, the evolved giant stars, and stars with 
a lot of activity. TurnbuU (2004) culled 131 such 
stars from a sample of 2300 within 30 pc concen- 
trating on "Sun-like" F, G, and K stars. A more 
extensive sample of stars suitable for the SETI 
program is given by TurnbuU & Tarter (2003) in- 
cluding approximately 18,000 candidates within 
2kpc. This "SETI sample" is important for its 
inclusion of stars at large distances that could be 
searched by very large space telescopes, and al- 
though they concentrate on solar-type stars, they 
include an incomplete sample of lower mass, a po- 
tentially restrictive limit as shown below. 



We examined each sample to find stars whose 
habitable zones can be resolved within the max- 
imum distance for telescopes of different diame- 
ters. The observational parameters for these cal- 
culations are: SN = 10, wavelength of 1 fiui, spec- 
tral resolution ^vjv = 0.01, detection efficiency 
rj ^ 1, observing time t = 24 hr. 

Table [1] lists the maximum distances and inner 
working angles derived under these assumptions 
along with the numbers of stars available for study 
from the SETI and 8pc samples. Table [2] shows 
how the numbers change if the inner working angle 
is assumed to be iX/dtei, typical of the best extant 
coronagraph designs, and allowing for an overall 
detection efficiency of 0.25. 

Figure [5] shows the cumulative number of stars 
whose habitable zones could be studied by sin- 
gle telescopes of different sizes for the SETI sam- 
ple. The bands of growth delineate the cumula- 
tive numbers using either the inner (smaller num- 
ber) or outer habitable zone radius. The numbers 
of candidates out to the maximum assumed dis- 
tances for the different telescope sizes are given in 
the figure as the average of these two criteria. 

For the SETI sample, the total number of ac- 
cessible stars increases by an order of magnitude 
when the telescope diameter doubles. This in- 
crease is more rapid than the d\^y rate predicted 
by the analytical calculations. A close examina- 
tion of the SETI sample shows that it deviates 
markedly from the PDMF in its distribution of 
stars, with a sharp peak in the number at about 
0.5 M0, a falling distribution toward lower masses 
and a flattening of the distribution between 1 and 
1.5 Mq,, completely unlike the local PDMF. This 
distribution over stellar masses also changes with 
distance, meaning the simple scaling laws derived 
for a constant PDMF will not describe the behav- 
ior of this sample. It is useful to note how these 
selection effects can distort the sample sizes avail- 
able for study, when additional criteria are used to 
select stars from the solar neighborhood that may 
be suitable for harboring life-bearing planets. 

These figures represent ideal observations, dis- 
counting operational inefficiencies, assuming the 
planet will be at maximum elongation at the time 
of an observation, ignoring additional background 
from the presence of interplanetary dust in the 
distant planetary system, and assuming perfect 
starlight elimination. If the interplanetary dust in 
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Fig. 5. — Cumulative number of stars from the 
SETI list vs. dtei using the ideal case of 77 = 1 and 
fe = 1, as in Table [T] Each telescope size shows 
sample growth in a band between two bounds cor- 
responding to 0IWA < Rin (lower bound) and 
OiwA < Rout (upper bound). The maximum de- 
tection distance from forms the right hand 
boundary of each sample. The numbers next to 
each diameter are the mean numbers of candidates 
between the two bands out to the maximum dis- 
tance, as appear in columns 4 and 5 of the table. 



distant planetary systems is typically as strong as 
it is locally, the effect will be to triple the amount 
of background light per observation on average in 
dH]): I{i^,z) 31{v,z), and the increased back- 
ground effectively decreases D^ax by a factor of 
(1/3)4 = 0.8, decreasing the sample size by about 
a factor of (1/3) V = 0.5. If 10"^ of the starlight 
leaks into the image of the planet the background 
increases by nearly an order of magnitude, further 
limiting the observing distance by a factor 0.56 
and decreasing the sample size by a factor of more 
than four . 

The inner working angle is a linear function of 
observing wavelength. We adopted A = 1 ^m as 
the minimum wavelength of interest to pick up 
spectral features, but it is obvious from Figure [T] 
that doubling or tripling this wavelength would be 
an enormous benefit to search for all the molec- 



ular signatures of life. Doubling the wavelength 
would decrease the sample size by another factor 
of about two for a coronagraph or an interferome- 
ter, where the inner working angle depends on the 
wavelength and the system size. 

Moreover, it is unlikely that a telescope whose 
inner working angle just resolves the inner habit- 
able zone will be adequate to observe exoplanets 
around any one star owing to the relatively small 
fraction of the projected orbit that will subtend 
large angul ar disp lacements for typical orbital in- 
clinations (jBrown 2 005: Br own ct al. 2 007). Re- 
alistically, OiwA will have to be much smaller than 
0IHZ to ensure a good chance of observing the 
planet in most cases. 

Each of these effects will decrease the sample 
sizes in the tables by factors of a few, and it is clear 
that the numbers in the tables might be more than 
an order of magnitude too large for the practical 
problem of life signatures. 

On the other hand, the range for terrestrial- 
mass planets thought to be suitable for life extends 
to 10 Mq, meaning the observed area of the planet 
in ^ will be about 10^^'^ « 5 times larger than 
7ri?|, , increasing the sample size for direct observa- 
tions by about a factor of 7 {N cx D^-^ cx R^-^ cx 
^2.5/3 ^ 102-5/3 « 7). This factor would boost 
the numbers in the tables to offset the decreases 
from other effects above. Furthermore, adopting 
24 hr as the limiting integration time could also be 
too conservative. For an integration time of 100 hr, 
the distance limit would increase by (35^)^/^, and 
the sample size by i^)'^'^^^ = 2.4. 

We can apply also ([48]) to these catalogs to find 
the number of stars where the atmospheric signa- 
tures can be detected through transit experiments. 
Table [3] gives the results. The I magnitudes were 
used to calculate the photon rates. The 8 pc sam- 
ple clearly limits the number of potential low lu- 
minosity stars for telescopes larger than a few me- 
ters. Using the analysis of §4.1) we can predict 
the sample of accessible stars (column 6) and the 
expected number from these samples where the or- 
bits align for transits (column 7). Evidently, the 
samples for transit observations will be orders of 
magnitude smaller than the samples for direct ob- 
servations for all telescopes. 

These uncertainties show why refinements to 
the parameters, for example using the greenhouse 
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effect to increase the size of the habitable zone 
slightly or trying to estimate the exo-zodiacal Ught 
background precisely, can complicate the analy- 
sis without bringing any greater insight into the 
likelihood of detecting life-bearing planets. More- 
over, the assumptions used here are optimistic 
about technological advances for observing exo- 
planets. There are no margins built in: we assume 
that all the starlight is rejected, there is no exo- 
zodiacal background, the detection efficiency is 
nearly ideal, and there is sufficient resolution when 
the inner working angle is equal to the size of the 
habitable zone. The numbers of candidate stars 
listed in Table [2] only take into account margins 
for 9iwA and detection efficiency, and we believe 
they may safely be taken as upper limits to the 
number of observable stars for a given telescope 
diameter, allowing for some margin in the uncer- 
tainties of the properties of exoplanets. Character- 
izing the atmospheres of Earth-like planets around 
other stars will be a challenging problem for the 
foreseeable future. 

6. Likelihood of Life Bearing Planets 

The number of life-bearing planets that can 
be detected in a survey will depend on the frac- 
tion of accessible candidate stars with an Earth- 
like planet in the habitable zones around the star 
as well as likelihood that the planets develop life 
quickly enough so that it starts to dominate the 
planet's atmospheric chemistry. This fraction is 
conventionally called ?7Earthl£l 

One estimate of the likelihood that a planet will 
be found in the habitable zone around a star can 
be derived from the fraction of stars with known 
exoplanets and the distribution of their orbital 
semi-major axes. There are presently 244 known 
exoplanet^, all much more massive than Earth, 
whose orbital parameters have been estimated. 
Figure [6] shows the distribution of semi-major axes 
of this sample. The range spans approximately 2.5 
orders of magnitude. In these logarithmic units, 
the habitable zone as defined by ^3]) is 0.27 wide, 
about 10% of the entire range. Because these ex- 



^There are various definitions of r^E^rth in the literature that 
may not require that the planet is in the habitable zone 
nor that it has developed life; the definition here is more 
demanding. 

^http://exoplanct.cu/catalog-all.php 



oplanets have been found in only 15% of the stars 
surveyed, the empirical chance of finding a mas- 
sive exoplanet in the habitable zone around a star 
is 1.5%. If this analysis also applies to terrestrial- 
mass planets, a sample of order 100 stars would be 
required to guarantee (Itr) at least one terrestrial 
planet in the habitable zone around a star. 



Distribution of 244 exoplanets 
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Fig. 6. — The distribution of the semi-major axes 
of 244 known exoplanets is shown along with the 
size of the habitable zone around the Sun. Al- 
though the location of the habitable zone for each 
candidate will differ from that of the Sun, the size 
will be the same in logarithmic units. 

This estimate could easily be too low for the 
problem at hand, however. There are many plan- 
ets in the Solar System, of which three, Venus, 
Earth, and Mars, might be capable of supporting 
life if they are within the habitable zone. The em- 
pirical chance of detecting massive exoplanets may 
be lower than the chance of finding a terrestrial 
planet, if the Solar System is typical. If planetary 
systems normally contain several planets of terres- 
trial mass, 0.5 - 10 Mq, the chance that one will be 
in the habitable zone rises considerably. The mi- 
crolensing detection of a 5.5 planet around an 
M dwarf (Beaufieu et al. 2006) suggests that ter- 
restrial planets may be more common around low 
mass stars than the empirical frequency of more 
massive planets. An optimist might say that ev- 
ery planetary system will have a planet within the 
habitable zone. The question of whether any could 
support life would depend on their mass; in the so- 
lar system, only about 1/3 of the known planets 
are in the right range. 

How likely is it that life has arisen and evolved 
to the stage where the atmospheric chemistry re- 
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fleets its presence? On Earth, life arose quickly 
and began to alter the atmosphere in ways that 
might have been detectab le from afar by an 
age of order 1 Gyr (jKaltenegger et al. 20071 ). It 
took more than 3.5 billion years to produce 
substantial amounts of oxygen via photosyn- 
thetic organisms to a state we would recognize 
today ( Kasting fc Catling 20031 and references 
therein). The oxygen would disappear from the 
atmosphere in few million years if life were to 
cease, and CO2 would dissolve in the oceans 
within a few thousand years, eliminating the most 
prominent atmospheric signature of organisms 
on Earth. It required a series of unlikely acci- 
dents for Earth to alter its atmosphere, suggest- 
ing that it may not occu r easily on other planets 
(IWard fc Brownlee 200GI) . 

Thus, an estimate made a priori of the fraction 
of stars with life-bearing planets depends entirely 
on how much faith we have that all the circum- 
stances are favorable, that is the number of ap- 
parent miracles we are willing to believe. If all 
stars have planetary systems (but only 15% have 
massive planets), and all planetary systems have 
at least one Earth-like planet within the habitable 
zone, and this planet always evolves life to domi- 
nate its atmospheric chemistry, then nearly 100% 
of the suitable stars will have planets indicative of 
life, i.e. ryEarth ~ 1- Even a 4m space telescope 
might be adequate to carry out the first observa- 
tions. 

On the other hand, if only 15% of stars typically 
have planetary systems of which the likelihood of 
an Earth-like planet in the terrestrial zone is only 
30% and the chance of evolving life is less than 
1, then 77Earth < 0.05, and it may be much less 
than this value if the sequence of events leading 
to life-signatures in the atmospheric spectra are 
more rare than common; many of the low mass 
stars making up the majority of the complete 8 pc 
sample are likely to be too young t o have evolved 
life as on Earth ijReid et al. 20071 ). A pessimist 
might conclude that jyEarth << 0.01. 

We believe that a pragmatic approach to the 
study of life-bearing exoplanets will require more 
than 100 candidate stars to yield at least one 
with the characteristics we seek, requiring a large 
(> 8 m) space telescope. A very large space tele- 
scope of order 16 m diameter would have thou- 
sands of candidate stars to study and, while tech- 



nically challenging to build, would also be an ex- 
cellent tool for examining exoplanets found by 
other means, such as spectra-photometry of the 
transits discussed in fjJl 

The range of uncertainty will narrow consider- 
ably when the results of the Kepler mission to find 
Earth-like planets around distant stars are known 
in a few years. However, Kepler is observing more 
than lO'^ distant (> 1 kpc) stars, very few of which 
will be close enough to look for atmospheric sig- 
natures. In the absence of very large space tele- 
scopes to find and study nearby stars, it will still 
leave open the question about about the likelihood 
of life outside of the Solar system. 

I am grateful to an anonymous referee who sug- 
gested several qualitative changes in the approach 
to this paper, vastly improving it. 1 also thank 
Robert Brown, Peter McCuUough, James Pringle, 
Marc Postman, Neill Reid, Kailash Sahu, David 
Soderblom, and Jeff Valenti for their advice and 
to Mike Hauser and Matt Mountain for encourage- 
ment to pursue this work. This research was sup- 
ported by NASA through its contract to AURA 
and the Space Telescope Science Institute and the 
University of California. 
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Table 1 

Sample sizes with ideal telescopes'" 



dtc\ 


dlWA 


^max 


n(SETI) 


n(8pc) 


m 


mas 


pc 


mz^ 


OHZ^ 


mz^ onz^ 


2 


123 


10 


2 


7 


4 8 


4 


62 


20 


26 


56 


9 16 


8 


31 


41 


300 


725 


19 29 


16 


15 


81 


3607 


6272 


30 55 



'^ejwji = 1.2^ @ 1 /im, = 1, < 24hr, ^ = 0.01 
""'OiwA < RpilHZ) or Rp{OHZ) 



Table 2 

Sample sizes with practical telescopes' 





^max 


n(SETI) 


n(8pc) 




(n(tot))'= 


m mas 


pc 


IHZ 


OHZ 


IHZ OHZ 


OHZ 


2 309 


6.8 








1 


1 


4 


4 155 


13.5 


3 


10 


3 


7 


35 


8 77 


27 


18 


78 


7 


14 


280 


16 39 


54 


241 


1092 


14 


22 


2240 




@ l^m, 


7] = 0.2 


5, t = 


24 hr, ^ = 


0.01 





^Calculated from ^7} 
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Table 3 

Sample sizes for transit detections'^ 





n(SETI) 


n(8p 






(^expect) 




-Rqhz 


^IHZ'^ 


^OHZ 


^IHZ 


Tp — 7© 


Tp = Tq 


2 





1 


3 


19 


8 





4 


3 


33 


22 


71 


62 


1 


8 


40 


256 


75 


87 


500 


5 


16 


284 


942 


87 


90 


4000 


40 



""Av/iy = 0.01, lOcr, 5yr, p = 5.5 gem ^, /i = 30amu, x = 4.2 

'^Numbers restricted by 8 pc sample limit 

'^Planet orbit at inner or outer HZ: Tp — 373 or 273 K 

'^Calculated from 

^Calculated from ([55]) 
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